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In this work the complete valence-band structure of the molybdenum dichalcogenides M0S2, 
MoSe2, and a-MoTe2 is presented and discussed in comparison. The valence bands have been studied 
using both angle- resolved photoelectron spectroscopy (ARPES) with synchrotron radiation, as well 
as, ab-initio band-structure calculations. The ARPES measurements have been carried out in the 
constant-final-state (CFS) mode. The results of the calculations show in general very good agreement 
with the experimentally determined valence-band structures allowing for a clear identification of the 
observed features. The dispersion of the valence bands as a function of the perpendicular component 
k± of the wave vector reveals a decreasing three-dimensional character from M0S2 to Qf-MoTe2 which 
is attributed to an increasing interlayer distance in the three compounds. The effect of this k± 
dispersion on the determination of the exact dispersion of the individual states as a function of fcy 
is discussed. By performing ARPES in the CFS mode the fe[| -component for off- normal emission 
spectra can be determined. The corresponding fcx-value is obtained from the symmetry of the 
spectra along the PA, KH, and ML line, respectively. 
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I. INTRODUCTION 

The molybdenum dichalcogenides belong to the large 
family of layered transition metal dichalcogenides whose 
crystal structure results from the stacking of sheets of 
hexagonally packed atoms. It consists of weakly coupled 
sandwich layers X-Mo-X in which a Mo-atom layer is en- 
closed within two chalcogen layers (X — S, Se, Te). All 
three M0X2 compounds crystallize in the 2Hb structure 
and M0S2, as well as, MoSe2 also crystallize in the 3R 
structure. In both of these polytype structures the coor- 
dination of the metal atoms is trigonal prismatic. In this 
work, 2IIb samples are investigated. For MoTe2 a phase 
transition from the semiconducting 2Hb-type a-MoTe2 to 
a metallic polytype, /3-MoTe2, with a monoclinic struc- 
ture having a distorted octahedral coordination is known 
Ma- The 2Hb-polytype, in the following assigned as 
MoTe2, is stable below « 815°C. 

The M0X2 compounds are indirect semiconductors 
with indirect (direct) band gaps of 1.29 (1.78) eV for 
M0S2, 1.10 (1.42) eV for MoSe2, and 1.00 (1.10) eV for 
MoTe2 0. Since the optical band gaps are matching 
well with the solar spectrum these materials are used for 
electrodes in high efficiency photo-electrochemical (PEC) 
cells J4J. In the molybdenum dichalcogenides the photo- 
transitions involve non-bonding d orbitals of Mo atoms. 
Therefore, these materials can be expected to resist hole- 
induced corrosion [g|. Despite of the major importance 
of the M0X2 compounds in solar cell production, to date 
no comparative band structure determination from ex- 



periment and theory has been reported. 

In general, it is expected that the gross features of 
the band structures of all three dichalcogenides are quite 
similar but, in addition, also some characteristic differ- 
ences are to be expected. Several band structure cal- 
culations were performed so far which were used to ex- 
plain experimental results from absorption measurements 
[^-|[ and early photoemission work ||10|-[l3|. From the 
literature about these materials no consistent picture 
emerges (see e.g. [ |l4| for a review, and references therein). 
Therefore, we address in this work the electronic band 
structure of all three binary dichalcogenides on equal 
footing by ARPES measurements and by 'state-of-the- 
art' ab-initio band structure calculations. In experi- 
ment, high-resolution photoelectron spectrometers (with 
respect to emission angle and energy) are used. In ad- 
dition, the constant-final-state (CFS) mode is apphed to 
obtain data which shows low k± dispersion. As to the 
theory, the electronic properties of all three molybdenum 
dichalcogenides are calculated within density functional 
theory (DFT) employing the local density approximation 
(LDA) [ll^ and spin-orbit interaction is included. 

The present photoemission data and ab-initio results 
constitute the first comprehensive picture of the elec- 
tronic valence-band structure of all three molybdenum 
dichalcogenides . 



II. CRYSTAL STRUCTURE 



The 2Hb-polytype of the crystal structure of the M0X2 
layer compounds is characterized by a stacking sequence 
/AbA BaB/ (A, B: chalcogen atom layers; a, b: Mo atom 
layers) p6[ . A unit cell and the corresponding first Bril- 
louin zone are displayed in Fig.^. fn this structure, a Mo 
atom layer is sandwiched between two chalcogen layers. 
The two X-Mo-X sandwich layers per unit cell are later- 
ally displaced with respect to each other so that the Mo 
atoms of the upper sandwich layer are directly above the 
chalcogen atoms of the lower sandwich layer in the unit 
cell and vice versa. The sandwich layers are coupled only 
by weak van der Waals forces giving rise to the quasi-two- 
dimensional character of the electronic structure. 

The lattice parameters of the direct lattice a, c, w and 
z are indicated in Fig.nl In this work, we use averages 
of the a and c values, as reported in the literature from 
different experiments P,n4,n6l . The thickness of a sand- 
wich layer is taken to be 2 z with z ~ 0.129 • c. Therefore, 
an interlayer distance w — c/2 — 2z — 0.242 • c follows. 
In Tab. | the used average experimental structure pa- 
rameters, the resulting c/a ratios, and the calculated re- 
ciprocal lengths of the most important high-symmetry 
lines (see Fig. |l| (b)) are given. In addition, the cova- 
lent radii of the chalcogen atoms and the ionic radii of 
the twofold negatively charged chalcogen ions are given. 
In the dichalcogenide crystals, the actual charge of the 
chalcogens is between the two limiting values of Oe and 
-2e as a simple population analysis of the calulated den- 
sity of states shows. It is obvious from Tab. || that the 
interlayer distance w, as well as, the lattice parameters 
a and c increase from M0S2 to MoTe2 mainly due to the 
increasing ionic radius of the chalcogen ions, while the 
c/a ratio remains nearly constant. 
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FIG. 1. Crystal structure (a) and the first Brillouin zone 
(b) of the 2Hb structure of the M0X2. 



TABLE I. Structural parameters of the 2Hb polytype of 
the molybdenum dichalcogenides for the direct and reciprocal 
lattice. The parameter z, i.e. half the thickness of a sandwich 
layer, is given by 2; = 0.129 • c. The size of the v. d. Waals 
gap w is calculated as c/2 — 2 2. Furthermore the size of the 
Brillouin zone is determined by the respective fc|| values. In 
addition, the covalent and ionic radii of the chalcogen atoms 
are given. The corresponding radii of molybdenum are Mo: 
1.40 A and Mo^+i 0.70 A, respectively. 





M0S2 


MoSe2 


MoTe2 


a [A] 


3.160 


3.299 


3.522 


c[A] 


12.294 


12.938 


13.968 


2. [A] 


3.172 


3.338 


3.604 


w [A] 


2.975 


3.131 


3.380 


c/a 


3.891 


3.922 


3.966 


FA [A-1 


0.255 


0.243 


0.225 


TM [A '] 


1.148 


1.099 


1.030 


FK [A-^] 


1.325 


1.270 


1.189 


MK [A '] 


0.662 


0.636 


0.594 


rcov [A] 


S: 1.04 


Se: 1.17 


Te: 1.37 


Vion [A] 


S^-: 1.84 


Se^-: 1.98 


Te^": 2.21 



III. EXPERIMENTAL APPROACH 
A. Sample preparation and experimental setup 

Naturally grown molybdenite (M0S2) was used for this 
work. Excellent samples with sizes up to 3 x 3 mm^ have 
been extracted from a larger piece. Single crystals of 
MoSe2 and a-MoTe2 were grown using the chemical va- 
por transport method (CVT) where bromine was used as 
the transport agent. 

Coarse orientation of the samples was achieved by 
LEED-patterns which were recorded under ultrahigh- 
vacuum conditions before cleavage. From the sharp- 
ness of the diffraction patterns it can be concluded that 
the single crystals are of high quality. Next, a small 
aluminum post was glued on top of the samples and 
they were inserted into the photoemission chamber. Af- 
ter preadjustment, the crystals were cleaved in situ by 
knocking off the aluminum posts. This was done at 
room temperature under ultrahigh-vacuum conditions 
(p< 2 X 10^ ^'^ mbar) to minimize surface degradation or 
contamination. Immediately after the cleavage a refer- 
ence spectrum was recorded to check the surface quality 
from time to time. Even after up to fourty hours of data 
acquisition at room temperature no significant changes 
were observed. This is due to the crystals cleaving at the 
van der Waals gap resulting in a rather inert surface. 

The measurements have been performed at the syn- 
chrotron radiation centers HASYLAB in Hamburg, Ger- 
many, and at BESSY-I in Berlin, Germany. The typical 



photon energy range of /ii/ = 9 ... 30 eV and a monochro- 
mator resolution of below 90 meV was used. The photo- 
electron spectrometers are hemispherical deflection an- 
alyzers mounted on two-axis goniometers. These are 
the WESPHOA-III station at the HONORMI beamhne 
(W3.2) at HASYLAB and the new photoemission sta- 
tion AR 65 at the 3m-NIMl and 2m-SEYA beamhnes of 
BESSY-I ||l^. For the measurements presented here, a 
monochromator resolution of 30-90 meV and an analyzer 
resolution of 50 meV were chosen, yielding an overall 
resolution of less than 100 meV. This is sufficient since 
all measurements were performed at room temperature. 
The maximum acceptance angle of the photoelectrons 
was« A^±0.8°. 



B. The CFS mode in photoemission 

Photoemission spectroscopy can be performed in sev- 
eral different modes to determine the occupied band- 
structure. Most common is the energy distribution curve 
(EDC) mode in which the photon energy is held con- 
stant and the detected kinetic energy is varied. With the 
emission angle i9 (with respect to the surface normal) 
and the detected kinetic energy Ekin of the photoelec- 
trons the wave-vector component parallel to the planes 
is determined by the dispersion relation: 
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(1) 



The fcj^-component is not conserved since the electrons 
travel through the crystal surface. If the final state band 
can be approximated by a free electron parabola starting 
at the inner potential Vq, the following equation can be 
deduced: 



\k±\ = j^/2m,{Ekin-cos^^ + Vo) 
n 



(2) 



Obviously, k± has to be larger than ^^/2meVo in order 
to be able to observe photoemitted electrons outside the 
sample. A key question for the analysis of the experi- 
mental data is whether this approximation can be used 
for the molybdenum dichalcogenides in the final state 
energy range used in this work. Recently Strocov et 
al. |lq,E9[ have claimed that this approximation might 
not be suitable for two-dimensional layered transition- 
metal dichalcogenides. We will come back to this point 
In addition, Strocov et al. 
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found for TiS2 
and VSe2 that the inner potential Vq is dependent on k for 
these materials. In the present study of the molybdenum 
dichalcogenides, Vq has been determined by LEED-I/V 
curves where intial electron energies between 50. ..250 eV 
have been used. For Vq we have determined 12.0 ± 2.3 
eV in M0S2, 14.5 ± 1.5 eV in MoSca and 16.1 ± 1 eV in 
MoTe2. An error of about 1.5 eV can be assumed for the 
determination of VJj by the LEED-I/V curve technique. 



Larger errors result from the statistical spread over sev- 
eral samples which have been investigated. For normal 
emission the critical points can be found accurately from 
the symmetry of the dispersion of the uppermost band 
along FAF. By this procedure the experimental band 
structure at the high symmetry points can directly be 
compared to the corresponding calculated band struc- 
ture. 

The constant-final-state mode has some advantages 
with respect to the EDC mode. For example, Eqs. y andg 
show that both fc-components depend on the kinetic en- 
ergy. Therefore, the emission features in an EDC spec- 
trum are affected by a dispersion in both directions. For 
a strongly dispersing band the corresponding peaks are 
asymmetric within the EDC spectra. This can be avoided 
in the CFS mode, where both fc-components (fcy and k±) 
are kept constant. A spectrum at a new value of e.g. fcy 

is adjusted by choosing a couple (Ekinjf?) for which k±_ 
remains constant. The range of these values is limited 
by the monochromator energy range so that not all high- 
symmetry directions are accessible. Another advantage 
of the CFS mode lies in the fact that the detected ki- 
netic energy at the analyzer is kept constant, too. Since 
the transmission function increases for lower Ekin values 
an EDC spectrum typically reveals an increased back- 
ground towards the low energy region, which cannot be 
determined exactly. In the CFS mode, only the charac- 
teristics of the monochromator (i. e. photon flux per pho- 
ton energy) influences the background of each spectrum. 
But this can be determined exactly, e.g., by measuring 
the current of a reference gold-mesh which is positioned 
in the photon beam. 



C. Data analysis technique 



Each spectrum measured in the present work has 
been modeled by a curve fitting program based on the 
Levenberg-Marquardt algorithm |2^]. The model func- 
tion is created by convoluting a sum of peaks plus Shirley- 
background pl| with the Gaussian shaped spectrometer 
function. Many emission features with larger binding 
energy reveal a Lorentzian line shape which is in agree- 
ment with the physical nature of the electronic states. A 
typical full width at half maximum of 200 . . . 400 meV 
indicates that the intrinsic line shape cannot be influ- 
enced by the spectrometer function since the latter was 
always chosen to be lower than 100 meV in width (i.e., 
analyzer plus monochromator resolution). Since in the 
CFS mode the whole spectrum belongs to one specific fc 
vector, the line shape cannot be influenced by dispersion 
effects which may occur in EDC spectra. 



IV. BAND-STRUCTURE CALCULATIONS 

The LDA-DFT calculations have been carried out em- 
ploying the norm-conserving, nonlocal pseudopotentials 
of Bachelet, Hamann, and Schliitcr [E2l. The exchange- 
correlation energy was taken into account using the 
Ceperley- Alder |23| form as parametrized by Perdew and 
Zunger |Q. As a basis to represent the wavefunctions, 
we use 120 Gaussian orbitals of s, p, d, and s* symmetry 
per sandwich layer and spin. At each atomic position, 40 
orbitals are localized. The decay constants of the Gaus- 
sians employed are {0.17, 0.45, 1.18, 2.00} for Mo, {0.17, 
0.45, 0.90, 1.70} for S, {0.17, 0.45, 0.90, 1.80} for Se, and 
{0.17, 0.43, 1.05, 2.60} for Te (in atomic units}. A real 
space mesh with linear spacings of about 0.2 A is used 
for the representation of the charge density and the po- 
tential. The spin-orbit interaction is treated in an on-site 
approximation, i.e., only integrals with the same location 
of the Gaussian orbitals and the spin-orbit potential are 
taken into account. Brillouin-zone integrations have been 
carried out using 12 special k-points in the irreducible 
part of the Brillouin zone. 

In our calculations we find, that a number of elec- 
tronic features depend sensitively on the lattice parame- 
ters. Therefore, we have used in the calculations reported 
in this work a value of z = 0.129 • c for half the diameter 
of the sandwich layers as obtained in experiment in or- 
der to arrive at a more meaningful comparison between 
theoretical and experimental results. Our calculations of 
z by total energy minimization (with a and c being held 
constant) yield 0.125 • c for M0S2, 0.128 • c for MoSe2, 
and 0.127 • c for MoTe2, i.e., values which are somewhat 
lower than the experimental values. 



V. RESULTS AND DISCUSSION 
A. Final states 

As mentioned above, the knowledge of the final state 
and the corresponding k± component is very impor- 
tant when the CFS mode is applied in photoemission. 
Eq. g assumes free-electron-like final-state bands which 
in general, as claimed by Strocov et al. [|l8|, might not 
be a good approximation for layered transition-metal 
dichalcogenides. As a result from their study of TiS2 
and VSe2, the authors argued that the complicated dis- 
persion of the final state bands affects the band mapping 
in a way that the k± component remains undetermined. 
In contrast, it can be expected for quasi- two-dimensional 
materials and especially for the M0X2 compounds that 
the E{k±) dependence is significantly smaller than the 
E{k\\) dependence. 

To determine experimentally the critical points at the 
zone edge a number of CFS spectra have been recorded. 
In particular, CFS spectra with fcii vectors along the FK 



direction at four different k± values corresponding to (0.0 
tt/c, -0.25 tt/c, 0.50 tt/c and 1.0 tt/c) were recorded and 
analyzed. The first of these four values corresponds to 
spectra taken along the FK and the last along the AH 
direction in the Brillouin zone. The experiments were 
carried out for M0S2 at the HONORMI beamlinc in HA- 
SYLAB (DESY, Hamburg) where a photon energy range 
oi hv = 10 ... 30 eV can be used. M0S2 was chosen since 
this material is expected to reveal the most pronounced 
three-dimensional character in the band structure. Nev- 
ertheless, the respective dispersion in most cases was 
found to be very small, as compared to the dispersion 
of the bands along the k\\ in-plane directions. Within the 
available photon energy range it was possible to record 
CFS series along the FK (AH) directions at four different 
k I values. 
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FIG. 2. The experimentally derived band structure of 
M0S2 along the FK and AH, as well as, two other parallel 
lines as indicated in the inset. 

Fig. presents the experimental band-structure fea- 
tures as resulting for E(fc||; fcj_ fixed) with the four dif- 
ferent fcj^ values noted above and indicated in the figure. 
For the sake of clarity only pronounced emission features 



in the spectra are shown in the figure. The filled cir- 
cles mark the FK direction with the k± value of the F- 
point. In general, the symbols of the other three series of 
measured features for the other three k± values lie very 
close to the circles. This is especially striking for the up- 
per bands near K/H. Only the two uppermost measured 
bands near F/A show a significant dependence on the k± 
value. Also the bands between -2.0 eV and -6 eV do not 
reveal a very pronounced E{k±) dependence. 




FIG. 3. Experimentally derived band structure (symbols) 
of M0S2 along the FK line in comparison with the respective 
calculated bulk band structure (solid lines) . 

At this point it is very revealing to compare some of 
the above discussed experimental results with results of 
our band-structure calculations. A respective compari- 
son for fcj_ equivalent to zero, i.e., along the FK direction 
is shown in Fig. 0. An amazingly good agreement be- 
tween theory and experiment for most of the bands is 
obvious in Fig. ^ lending further support to the accurate 
experimental determination of the high symmetry line. 
Most of the Mo Ad- and S 3p-derived bands between 
eV and -5 eV are found to be in very good agreement. 
Only the S 3p-derived band around -6 eV shows signif- 
icant deviations. One possible origin of this behaviour 
in experiment could be the large width in the underlying 
data, especially at F, so that the assignment of the mea- 
sured structures to a well-defined binding energy is diffi- 
cult. Another origin of the seeming disagreement could 
be related to the symmetry dependence of the observed 



spectral features. We come back to this point in Sec- 
tion V B. The band in question originates from sulfur 
3p valence states and its dispersion is determined by the 
interaction of these states across the van der Waals gap. 
To obtain further evidence that the K and the M point 
at the zone edge can be accurately determined, CFS 
series along the KH and ML high-symmetry lines were 
recorded. The respective k± values have been calculated 
by using Eq. even if this approximation is quite coarse. 
This is sufficient since the determination of the critical 
points has been done experimentally. In this case, MoSe2 
was chosen as the example since it shows a distinct spin- 
orbit splitting along the KH line in addition to a clear 
three-dimensional band character. After analyzing the 
spectra by the curve fitting approach discussed above, a 
weak E{k±) dependence is observed in the data. 
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FIG. 4. Experimental valence-band structure (symbols) of 
MoSe2 along the symmetry lines KH and ML in comparison 
with the calculated bands (solid lines). 

This is shown in Fig. ^ where the experimental data 
are compared to the corresponding calculated bands. 
The energy axis is aligned to the valence-band maximum 
which results from further investigations presented in the 
next sections. First we note that there is a reasonable 
general correspondence between the number and energy 
position of measured and calculated bands. Even the dis- 
persion of the theoretical and experimental results show 
good agreement in cases. The absolute binding energy for 
the topmost emission features, however, are slightly lower 
than the energies of the respective calculated bands. One 
can therefore conclude from the dispersion and energy 
position of the topmost bands and from the bands in 
the energy region around -2.5 eV that the K and the M 
points can be determined. By having the three high sym- 
metry points, namey F, K, and M, one is able to scan all 
high-symmetry lines. 

Despite a good correspondence of the states with 
smaller binding energies up to -3.0 eV the states between 
-3.5 eV and -4.5 eV along the KH symmetry line show a 
fairly poor agreement with theory. In particular, the cal- 



culated bands show a significant k± dispersion while ex- 
periment observes nearly dispersionless states. We have 
observed the same behaviour already in Fig. g, where two 
bulk bands at about -6 eV and -7 eV were found while ex- 
periment observed only one very narrow band. We come 
back to this point further below. 

For the ML direction a direct comparison of the bands 
in the latter energy region is not possible, since the re- 
spective bands could not be observed in experiment be- 
cause of the limited monochromator energy range. 



B. Valence-band structure along FA 

Normal emission spectra for all three materials were 
recorded at various beamlines at HASYLAB or BESSY- 
I synchrotron sources, respectively, in the energy range 
of 13. . .30 eV. The energy axis is referred to the valence- 
band maximum which was determined in connection with 
the complete valence-band structure as described in the 
following sections. To determine the band structure from 
the spectra series all spectra were modeled after the re- 
moval of a small Shirley background pi] by a sum of 
Gaussian and Lorentzian profiles. The k± value within 
the Brillouin zone to which each peak has to be attributed 
was calculated from Eq. 0. As pointed out above, in this 
limited energy range of about 10. . .30 eV this equation 
can be used allowing for a sufficiently precise calculation 
of the k± value for each energetic peak position. For 
larger energies the deviation of the final states from free- 
electron-like behavior must be taken into account. Here 
target current spectroscopy might be applied to first de- 
termine the coupling bands and then measure all pho- 
toemission spectra in the CFS mode with the final state 
belonging to the specific coupling band ^9 25,2^]. 

Before addressing the observed band structures along 
FA in detail we note that only each second F point is ob- 
served in experiment. According to Ref. u% this is due 
to the fact that in non-symmorphic crystal structures the 
symmetry of the final state is not always unique as it is 
for symmorphic structures. The FA direction A con- 
sists of the points k = 27r/c(0,0, i5) with < S < ^. A 
represents a group with a number of operations which 
are listed for instance in the above reference. For exam- 
ple a plane wave of the form exp[z27r/c(0, 0, <5)(x, y, z)] is 
transformed by any operation from A to an irreducible 
representation and this plane wave belongs to e.g. Ai. In 
contrast a plane wave exp[z27r/c(0, 0, 6+l){x, y, z)] which 
is streched by a reciprocal lattice vector 27r/c(0,0, 1) be- 
longs to A2. Obviously the symmetry depends on the 
length of the wave vector. For a given final state symme- 
try this leads to the observed missing of each second F 
point since each second initial state reveals an odd sym- 
metry. This fact has to be taken into account in the 
comparison to the band structure calculations. 



our band structure calculations. In one case (Fig. ||) we 
have included the density of states of the bands along 
the FA direction, as well. We first address the experi- 
mental results. In general, for all three materials one can 
easily distinguish at least two dispersive structures and 
a number of nearly non-dispersive bands. The upper- 
most valence band, as well as, the band with the largest 
binding energy reveal a significant fc_L-dependence. The 
width of the uppermost band between F and A decreases 
from M0S2 to MoTe2, i. e., the bands become more two- 
dimensional and less fcj_-dependent. This effect is not 
so strong for the band with the largest binding energy 
but also noticeable. From this result one may conclude 
that the latter might be predominantly of chalcogen pz 
orbital character, which is relatively strongly influenced 
by the following sandwich layer in the crystal. Indeed, 
our band structure calculations exhibit a contribution 
of the chalcogen pz orbitals from 66% to 60% for the 
three chalcogenides to this state. In contrast, the influ- 
ence of the next sandwich layer on the predominantly Mo 
4dj,2-derived states (61% to 68% metal d^2 character) at 
the valence-band maximum decreases rapidly with the 
increasing size of the unit cell. 
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The experimental band structures along FAF are 
shown in Figs, g, ^, and together with the results of 



FIG. 5. Experimental valence-band structure (symbols) of 
M0S2 along the FA high-symmetry line in comparison with 
the calculated band structure (solid lines). Open symbols de- 
note weak measured structures. The error bars are due to 
the limited experimental energy and angular resolution and 
due to uncertainties from the curve fitting. In addition, the 
density of states for the bands along FA is plotted. The con- 
tributions of particular molybdenum d-states is indicated by 
the shaded regions. 
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FIG. 6. Experimental valence-band structure (symbols) of 
MoSe2 along the FA high-symmetry line in comparison with 
the calculated band structure (solid lines). Open symbols de- 
note weak structures. The error bars are due to the limited 
experimental energy and angular resolution and due to un- 
certainties from the curve fitting. 
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FIG. 7. Experimental valence-band structure (symbols) of 
MoTe2 along the FA high-symmetry line in comparison with 
the calculated band structure (solid lines). Open symbols de- 
note weak structures. The error bars are due to the limited 
experimental energy and angular resolution and due to un- 
certainties from the curve fitting. 



It should be noted that in all compounds the topmost 
band at the second F-point (for which it can be observed) 
lies at larger binding energies in comparison to the first 
one. This may be attributed to two effects. First, the an- 
alyzer has a fixed maximum acceptance angle 'dmax ■ The 
fc-resolution Afc depends on this acceptance angle and 
on the kinetic energy (see Eq. nl). Assuming ii}„iax — 0.8° 
then for the first F-point Ak « 0.020 [A^i] while for the 
second one it is Ak « 0.036 [A^^]. Since the F-point is 
the point with the highest symmetry the band disperses 
in all fc-directions away from F to larger binding ener- 
gies. For a lower fc-resolution the band position at all k 
values within this range must be integrated. This yields 
an emission feature with a maximum at a larger binding 
energy. The second effect might be the reduced informa- 
tion depth at higher photon energies. If the mean free 
path becomes lower than one unit cell, the k± depen- 
dence will smear out and the spectra will reflect a more 
one-dimensional density of states |Q . 

Below the two dispersive bands a number of nearly 
perfect dispersionless bands are observed. The density of 
states resolved with respect to Mo M and S 3p-character 
in Fig. |5| clearly reveals the physical origin of some of 
these features. From the results of our band-structure 
calculations it is most obvious, that these bands are two- 
dimensional and can be attributed to Mo Ad states (see 
the calculated density of states in Fig. g for that matter). 
In the calculations, they are located at about - 1.8 eV and 

- 3.1 eV and the bands are separated in both cases due 
to spin-orbit interaction. For M0S2 this splitting is too 
small to be observed in experiment. In contrast, for the 
other two compounds the photoemission data clearly ex- 
hibit the splitting of the dispersionless bands at about 

- 1.5 eV and - 1.1 eV, respectively. In MoSe2 this ener- 
getic separation is about 300 meV and for MoTe2 it is 
about 500 meV. So, the energetic separation of these two 
spin-orbit splitted band groups, as well as, their binding 
energy decrease towards the tellurium compound in very 
good agreement between theory and experiment. The 
density of states in the other energy regions (see Fig. ||) 
is comparatively small, so that respective band-structure 
features are more complicated to measure. 

Thus a comparison with the calculated band structures 
in general shows very good agreement, indeed, in spite of 
the fact that some experimentally observed emission fea- 
tures do not result from the calculations. In all three 
compounds such emissions are found above the lowest 
binding energy state at about -4.5 eV. Typically this is 
only a very weak structure but it is also directly visi- 
ble in the spectra of e.g. M0S2 in the energy range of 
21. . .23 eV. Additional states may also arise from transi- 
tions of an initial state to a neighboring final state with 
different binding energy. Since an energy distribution 
curve provides only combined information of the initial, 
as well as, the final state it cannot be distinguished be- 
tween such effects. Of course, one would expect that only 
one coupling band is suitable for a given initial state but 



as it has been stated by Strocov et al. ||l8[ , in some cases 
an additional final state branch may exist. The emis- 
sion intensity is then expected to be lower as from the 
main branch and this can be supported by the observed 
features in the present spectra. 

Additional valence bands between the non-dispersive 
bands and the topmost valence band are experimentally 
observed but not supported by theory. A number of the 
respective spectral features are quite intense and might 
not be explained by different coupling branches. In M0S2 
at least two strong emission features appear at about - 
1.3 eV and -0.8 eV which are quite non-dispersive for 
lower photon energies. For higher photon energies, i. e., 
beyond the second F point the band near -0.8 eV begins 
to disperse to smaller binding energies while the band 
which is found energetically below the former vanishes. 

Finally, we note that all three valence band structures 
in Figs, ra, H and R show valence bands within the energy 
region from - 5 eV to - 7 eV with significant fc_L disper- 
sion. This dispersion shows the symmetry properties, as 
discussed above, in that only one of the two bands be- 
tween consecutive F-points is seen in experiment. The 



agreement between the respective theoretical and exper- 
imental bands is particularly close for MoTe2 as can be 
seen in Fig. ^ By the same token, it becomes most ob- 
vious, why experiment observes in this energy region for 
one particular F point cither the upper or the lower band 
but not both. Since the experimental data plotted in 
Fig. were recorded for one particular k± at the F-point 
only one band could be observed. 



C. Valence- band structure along the in-plane 
directions 

After the discussion of the valence-band structure in 
the FA direction we now focus on the directions parallel 
to the layers. We discuss the band structures along the 
high-symmetry lines FK, FM, and MK of the Brillouin 
zone. The FK direction is the longest high-symmetry 
line. In the case of the molybdenum dichalcogenides this 
direction is most important since the energetic position 
of the topmost emission features at the F and the K point 
determine a number of electrical and optical properties. 




MoSe 




FIG. 8. ARPES spectra recorded in the CFS mode along the high-symmetry line FK of the three molybdenum dichalcogenide 
compounds. The size of fcy at the K-point is given in Tab. |. 



Fig. g shows CFS spectra for the three molybdenum 
dichalcogenides recorded along the FK line. The respec- 
tive final state was chosen such that k± remains in the 
FKM plane within the whole series of spectra. Here the 
experimental determination of the F, K, and M point has 
been taken into account. 

A number of features is easily recognized in the spectra 
of Fig. g. The three series reveal a close similarity of the 
E{k) behavior of most emission features. Starting at F 
(i. e., fc|| = 0.0 A~^) to the K point the states near -3.2 eV 
and the lowest binding energy disperse towards larger 
binding energies. This is different from the two topmost 
emission features, which separate from each other after 
nearly halfway towards K and move with nearly constant 
width to K. 

From the spectra in Fig. ^ it can be deduced that 
for M0S2 the topmost emission feature near F reveals 
the weakest intensity compared to the other compounds. 
Since this state is predominantly derived from the Mo 
4:dz2 orbitals (see previous section) high spectral inten- 
sity for this state is expected. It should be decreased 
for a lower lattice parameter a due to a shielding by the 
chalcogen p orbitals. Indeed this is observed. Further- 
more, the topmost emission feature reveals a decreasing 
dispersion from M0S2 to M0TC2. 

The spectra for the FK and FM directions are largely 
similar. The spectra along the MK direction show rela- 
tively small dispersions, as compared to those along FK 
and FM. Only the features with the smallest binding en- 
ergy show appreciable dispersion. 
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FIG. 9. Double emission feature for the three dichalco- 
genides at the K point. Filled circles denote the experimental 
data while the solid lines represents the fit curve as obtained 
from the sum of two Gaussians (dashed lines) (X, Y) and an 
additional weak emission feature labelled '"add."'. 



D. The spin-orbit splitting at K 



Before discussing the complete valence-band structure 
of the three molybdenum dichalcogenides, we first briefly 
address the effects of spin-orbit coupling on the spectra 
and the band structures. A distinct splitting of the emis- 
sion features with lowest binding energy at the K point 
can be observed in Fig. |[ It is most obvious for MoTe2 
and is only visible as a weak shoulder in the M0S2 com- 
pound. Fig. H shows the relevant part of the spectra for 
each compound where the double emission feature is as- 
signed with X and Y, respectively, in line with previous 
work on MoTe2 [£9[ . A broad additional emission feature 
denoted as 'add.' is taken into account which is present 
in all three spectra just below peak X. 



The size of the splitting of these bands at the K-point 
in the experimental spectra and in the theoretical results 
is listed in Tab. ||. From M0S2 to MoTe2, the splitting 
increases with increasing Mo d character of the bands. 
Furthermore, the energetic position of the bands move 
towards the valence-band maximum at the F-point (cf. 
Fig. |). 

The general trend in the observed energetic differences 
whithin the row of the three dichalcogenides is observed 
in the theoretical results, as well. The experimentally 
observed splittings of the emission features X and Y at 
the K-point, however, are about 100 meV smaller than 
the respective calculated splittings. 

Another point concerning the splitting at K is con- 
nected with the interpretation of two structures in opti- 
cal measurements. Optical absorption spectra of all three 
M0X2 compounds display excitonic features character- 
istic of the group VIA transition metal dichalcogenides 
P,^ . As claimed by Coehoorn et al. | ll| , who observed 
a splitting of 0.21 eV for MoSe2 at the K point, these 



states give rise to the so-called A and B excitons. The 
splitting found in Ref. |11|| for MoSe2 matches well with 
the observed e nerg etic difference of the A and B excitons 
(see, e.g., Ref. [31]) and therefore the authors concluded 
that probably the exciton peaks arise from the two pos- 
sible transitions at the K point. 

TABLE II. Energetic separation of the emission features 
X and Y at the K point in comparison with the calculated 
splitting. The right column gives the energetic distance be- 
tween the A and B excitons following from optical absorption 
measurements la. 



these bands should be correct. The comparison with the 
results of inverse photoemission investigations on M0S2 
p5| , p6[ confirms this finding. Also studies on the simi- 
lar material WSe2 show good qualitative agreement [g^ . 
Thus it can be concluded that the minimum of the con- 
duction bands is located nearly halfway between F and 
K, indeed. Furthermore, the valence-band maximum is 
found at F in all three compounds. This is confirmed for 
all three materials by the experimental data. From the 
present band structure calculations it follows that the en- 
ergetically lowest indirect transition is from F to 
and the lowest direct transition occurs at K. 



M0S2 

MoSea 

MoTe2 



Splitting at K 
Exp. ARPES data Calculation 



161 ± 10 meV 
175 ± 10 meV 
238 ± 10 meV 



258 meV 
294 meV 
331 meV 



A, B excitons 
Exp. opt. data 



180 meV 
250 meV 
380 meV 



Additionly, Finteis et al. found the valence-band max- 
imum of WSe2 to be located at the K point in their ex- 
periments |£^,^ . The latter authors also concluded that 
the double emission feature at the K-point is the signa- 
ture of the two initial states of the excitons A and B. 
To attribute the double emission feature at the K point 
to the initial states of the excitons A and B is at vari- 
ance, however, with several investigations using optical 
absorption measurements pJtIJS^I leading to the result 
that appropriate transitions for the A and B excitons 
should occur at the F point. However, according to the 
present experimental data the peak distances as listed in 
Tab. O for the double feature X, Y for all three materi- 
als are too small compared to the energetic distance of 
the excitonic lines measured by e.g. Beal et al. M. Even 
though a shift between the initial states and the observed 
energy distance of the excitons is possible [Q the mag- 
nitude of the observed shift cannot be explained within 
the theory [Q. It must be concluded accordingly that 
the double emission feature X, Y at the K point does not 
correspond to initial states of the excitons A and B [g9| . 
In line with the optical absorption measurements cited 
above and from the location of the VBM at the F point 
it is suggested that the topmost occupied state at F is 
the initial state for the A and B excitons. 



E. Complete valence-band structures 



In Figs. 
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nu and n3 th e ex perimentally observed 
spectral features (see section |VC] ) of each molybdenum 
dichalcogenide compound are plotted as a band structure 
along the high-symmetry lines MFKM in direct compar- 
ison with the results of our band-structure calculations. 
A projected band structure has not been used since the 
respective critical points have been determined experi- 
mentally. 

The lower unoccupied bands are also shown. Except 
for the gap problem within the LDA the dispersion of 



iFK 



A close look at the valence-band structure of all three 
materials reveals a very good general agreement between 
measured and calculated bands. The total valence-band 
width is found to decrease from M0S2 over MoSe2 to 
MoTe2 both in experiment and theory. 

The dispersion of the uppermost bands away from F to- 
wards K and M reveal also quantitative matching. How- 
ever, the energetic positions of the calculated bands in 
the vicinity of the high-symmetry point K differ with a 
maximum difference of 200 meV in the case of MoSe2. 
Also the calculated splitting of the bands at K (already 
discussed above) is about 100 mcV larger than observed 
in experiment. The general trend in decreasing binding 
energy of the top of the valence bands at M, as observed 
in experiment, results from theory, as well. Only the 
absolute energetic location of the topmost state at M is 
about 200 meV too low as compared to the photoemis- 
sion results. It must be emphasized that such differences 
between experiment and theory in the range of a few 
100 meV are close to the theoretical uncertainty that can 
be achieved nowadays. 

In the vicinity of about -1 eV a dispersionless band is 
observed in the FA direction (see Figs. g| ^, and0) which 
is also dispersionless near the center of the Brillouin zone 
for off- normal emission (see Figs. HG, 11, and [l2|) [p8[. 



In the FKM-plane, this band is found experimentally in 
M0S2 to be «100 meV below a rather dispersionless band 
in the calculated band structure and about 100 meV 
above the respective band for MoSe2 and MoTe2. On 
the other hand, along FA only in the case of MoTe2 this 
band results from the theory. 

In band structure calculations for the (OOOl)-surface of 
M0X2 (X=S, Se, Te), which we have carried out using a 
slab of 9 sandwich layers (method as described in |p7|), 
we have found states located at the first few sandwich 
layers of the surface in this energy range. The respective 
bands exhibit only a small dispersion around the center of 
the surface Brillouin zone in agreement with the experi- 
mentally observed bands. Therefore, one might attribute 
these bands in the photoemission data to surface states. 
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FIG. 10. Experimental valence-band structure (symbols) 
of M0S2 along the in-plane directions MF, FK, and KM in 
comparison with the calculated band structure (solid lines). 
Open symbols denote weak structures. The error bars are due 
to the limited experimental energy and angular resolution and 
due to uncertainties from the curve fitting. 




FIG. 12. Experimental valence-band structure (symbols) 
of MoTe2 along the in-plane directions MF, FK, and KM in 
comparison with the calculated band structure (solid lines). 
Open symbols denote weak structures. The error bars are due 
to the limited experimental energy and angular resolution and 
due to uncertainties from the curve fitting. 



VI. CONCLUSIONS 




FIG. 11. Experimental valence-band structure (symbols) 
of MoSe2 along the in-plane directions MF, FK, and KM in 
comparison with the calculated band structure (solid lines). 
Open symbols denote weak structures. The error bars are due 
to the limited experimental energy and angular resolution and 
due to uncertainties from the curve fitting. 



In our investigations, we have determined the complete 
valence-band structure of M0S2, MoSe2, and a-MoTe2 
on equal footing by experiment and theory. It was found 
that the valence bands resulting from angle-resolved pho- 
toelectron spectroscopy in the constant-final-state mode 
and from ab-initio calculations are in very good overall 
agreement. Our results yield the first comprehensive pic- 
ture of the electronic bulk valence-band structure of all 
three molybdenum dichalcogenides. 
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